Non-adiabatic preparation of spin crystals with ultracold polar molecules 
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We study the growth dynamics of ordered structures of strongly interacting polar molecules in op- 
tical lattices. Using dipole blockade of microwave excitations, we map the system onto an interacting 
spin-1/2 model possessing ground states with crystalline order, and describe a way to prepare these 
states by non-adiabatically driving the transitions between molecular rotational levels. The pro- 
posed technique bypasses the need to cross a phase transition and allows for the creation of ordered 
domains of considerably larger size compared to approaches relying on adiabatic preparation. 

PACS numbers: 67.85.-d, 34.20.Gj, 42.50.Dv, 03.75.Hh 



Long-range dipolar interactions enable the creation of 
novel states of matter with ultracold quantum gases [T] . 
A prominent example is the predicted formation of dipo- 
lar crystals with tunable interaction parameters [3S]- 
However, the preparation of such strongly interacting 
phases starting from a weakly interacting quantum gas 
is known to be very challenging, as it involves crossing 
a phase transition, where the energy gap vanishes in the 
thermodynamic limit [5|. In this Letter we show that 
for ultracold polar molecules in an optical lattice, this 
problem can be overcome by non-adiabatic driving of ro- 
tational transitions. 

Our approach builds on recent experimental advances 
in the coherent creation and control of ultracold po- 
lar alkali metal dimers, such as KRb [6] and LiCs [7]. 
The experiments have demonstrated that ultracold po- 
lar molecules can now be produced in the rovibrational 
ground state, transferred to any hyperfine sublevel [S], 
and trapped in a periodic potential of an optical lat- 
tice [5] . Ultracold molecules trapped on an optical lattice 
have previously been proposed as promising candidates 
for quantum computation [lOj and quantum simulation 
of spin-lattice models [J [TTJ [T^] . Here, we employ dipole 
blockade of microwave excitations in the context of an 
effective spin-1/2 model with ultracold molecules on an 
optical lattice. The ground state phase diagram of such 
model is dominated by a series of commensurate phases, 
in which one of the spin states exhibits crystalline or- 
der. We show that these phases can be prepared by a 
short sequence of microwave pulses that nucleate the or- 
dered domain, followed by continuous microwave driving 
that propagates the domain boundary. We provide an 
effective model that describes the dynamics under the 
continuous driving and demonstrate that it leads to an 
efficient growth of the ordered domains. Finally, we an- 
alyze the imperfections and the required experimental 
parameters, and demonstrate that structures consisting 
of > 1000 spins can be grown. 

To explore these ideas, we consider a one-dimensional 
array of ^E molecules with a dipole moment d prepared 



in a Mott insulator state in an optical lattice with the 
period a, tilted by the angle with respect to the lab- 
oratory z axis, as shown in Fig. fl] (a). The lattice is 
filled with one molecule per site and the trapping poten- 
tial is strong enough to prevent tunneling between the 
sites and molecular reactions [9 . The rotational states 
of a ^E diatomic molecule, | J, M), are labeled by the an- 
gular momentum J and its projection on the laboratory 
z-axis, M . The energy of the rotational states is given 
by i?rot — BJ{J+1), where B is the rotational constant. 
We assume that the molecules are initially prepared in 
the rotational ground J = state. 

The states \J, M) can be coupled by a microwave field, 
producing linear combinations of rotational states with 
different parity. In particular, the rotational states |0, 0) 
and 1 1 , 0) can be coupled by near-resonant linearly po- 
larized microwave field with Rabi frequency O,, and de- 
tuning Ac to produce two field-dressed states separated 
by ~ ilc- If the microwave field is applied adiabatically. 
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FIG. 1. (a) ID array of polar molecules on an optical lattice 
tilted by the angle 6 with respect to the laboratory z-axis. (b) 
Lowest rotational states of a polar ^E molecule. A strong mi- 
crowave field Q.C couples the states |0,0) and |1,0), providing 
them with permanent dipole moments in the rotating frame; 
a weak field ^ <^Q.c drives the two-photon transition between 
one of the dressed states \g) and the state |e) = |3,0), with 
the single-photon frequency far detuned from the J — 2 state. 
All fields are linearly polarized along z. 



the molecules must all populate the same field-dressed 
state. We choose this state, \g) = a|0, 0) -l-6|l,0), as 
the ground, 'spin-down', state. In the rotating frame, 
the state \g) has a permanent dipole moment dg — 
{-J2ah/-j2,)d, leading to the dipole-dipole interaction be- 
tween the molecules, Vdd('') — dPg/r^i^ — 3cos^ 9), whose 
magnitude and sign are tunable by changing Ac/Slc and 
9 [13j. We assume the coupHng field to satisfy the con- 
dition Vdd(a) ^ r^c <C IB that ensures that the dipole- 
dipole interaction does not mix \g) with the neighboring 
field-dressed state nor with the dark |1,1) state, which 
for typical experimentally realizible systems (molecules 
with dipole moment 1-5 Debye and a ~ 250 - 500 nm) 
requires fie ~ 1 — 100 MHz. As the effective spin-up state 
we choose a level without a permanent dipole moment, 
|e) = |3,0). The two pseudo-spin states, \g) and |e), are 
connected by a weak two-photon transition with Rabi 
frequency fJ ^ f^c, with the single photon resonance far 
detuned from the |2,0) state, which remains unaffected 
by the microwave field. 

Using the rotating wave approximation, we can write 
the Hamiltonian for an ensemble of molecules on an op- 
tical lattice as 

H = -r.A^ |e,)(e,| + ^ 5^(|e,)(5.| + k)(e.|) 
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where Y = Vdd(a) and A = a; — Wge is the detuning of 
the two-photon field from the g — e resonance. Eq. (fTl) 
can be expressed via spin- 1/2 operators Sa as 



H = h^Sf ^h^S^ +VY^ 
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(2) 



where h^ = hn, h^ = -{HA + C(3)V"), and ((3) = 
J2T=i l/fc^ ~ 1-202 is Riemann's zeta function [11]. The 
magnitude and sign of the parameters h^/hx and V/hx 
can be tuned by changing A, 9, and Ac/^c- Note that 
the two-photon transition driven by fl is key to realiz- 
ing the dipole blockade [15] with molecular rotational 
levels: the pseudo-spin states \g) and |e) have a differ- 
ence in angular momentum of A J > 2, and therefore are 
not mixed by the dipole-dipole interaction. This leads 
to strong interactions between molecules in the \g) state, 
while eliminating the interactions between molecules in 
the \g) and |e) states (the "flip-flop" terms). 

The thermodynamic properties of the Hamiltonian ([2]) 
have been studied before in the context of Rydberg 
atoms |16H18j , and similar implementations based on po- 
lar molecules have also been discussed previously [19]. 
The main advantage of using polar molecules is that the 
dynamics occurs within the manifold of low-energy ro- 
tational states, whose lifetimes (> 1 s) are much longer 
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FIG. 2. Schematics of the ground-state phase diagram corre- 
sponding to the Hamiltonian H. The commensurate crystal 
phases with indicated fillings (orange) are surrounded by the 
floating solid (FS) phase (green). 



than the lifetimes of Rydberg states (~ /zs). Therefore, 
polar molecules offer the possibility of creating ordered 
structures with much larger size. 

The main features of the phase diagram of the system 
described by Eq. ^ are shown in Fig. |2] and can be best 
understood in terms of the flUing with which the minority 
spin component occurs. When V is positive and domi- 
nates over both h^ and h^, the system is in an antiferro- 
magnetic (AFM) phase, corresponding to half-fllling, i.e. 
the up and down spins occur in an alternating order. In 
the absence of a longitudinal field, h^ = 0, with V be- 
ing large and negative, all the spins are pointing in the 
same direction, which corresponds to the ferromagnetic 
(FM) ordering. There, any finite value of h^ breaks the 
Z2 symmetry of the model leading to the breakdown of 
the FM phase [20]. Note that the long-range interaction 
leads to the opposite shifts of the transition points on 
the hz = line. In the FM phase, the long-range terms 
lead to renormalization of the nearest neighbor coupling 
constant, while for the AFM interaction they result in 
frustration accelerating the melting of the phase. From 
both sides of the AFM phase there is a complete devil's 
staircase of crystalline configurations with different com- 
mensurate lattice spacings [TT] [T51 [H] . In general, the 
commensurate phases melt in two steps: first into the 
floating solid (FS) phase, and then into the paramagnet 
(PM) . The FS phase is characterized by a finite density of 
Bose condensed dislocation defects of the commensurate 
crystal and exhibits a gapless excitation spectrum [17| . 
At the tip of the lobe at filling 1/3 and 1/4, however, the 
system belongs to the Potts universality class, which re- 
sults in a direct transition from the commensurate phase 
tothePM HSj. 

Adiabatic preparation of the ordered phases, as e.g. 



proposed in Ref . [12] , is challenging due to the vanishing 
energy gap between phases on the phase boundary. For 
a system of N particles the gap closes as 1/A^ for the 
Ising transition between the PM and the AFM phases 
[22] and for the transition between the PM and the FS 
phases [53], and as 1/N'^ for the transition between the 
FS phase and the commensurate crystals due to its re- 
lation to free fermions [5l|T7]. Within the Landau-Zener 
approximation (equivalent to the Kibble-Zurek formalism 
for continuous phase transitions) for the adiabatic cross- 
ing of the phase boundary [22], the average size of the 
ordered domains is proportional to the square root of the 
sweep duration. While, in some cases, this scaling can be 
improved using non- linear sweeps [231 I24j , it is highly de- 
sirable to have a more robust protocol. Here, we present 
a non-adiabatic method for creating crystalline order in 
the system, with the growth of the domains scaling much 
better than within the adiabatic approach. 

In the following, we use the dipole blockade of 
microwave excitations to non-adiabatically construct 
strongly interacting ground state phases of the Ising 
model, eq. (l2]), as schematically shown in Fig. [SJ We 
consider the case of ^ > 0, corresponding to 9 — tt/2, 
and exemplify the technique by caclulations for ^Li^'^'^Cs 
molecules (d = 5.520 Debye, B = 5.636 GHz [2F]) in an 
optical lattice with a = 266 nm. 

The procedure to grow large ordered domains is 
schematically illustrated in Fig. [3] In the absence of the 
spin-flipping field f2, all the molecules are initialized in 
the ground (spin-down) state \g). The energy needed 
to resonantly flip the first spin in the lattice is given by 
El = —2V (,{?>), and we use it to renormalize the detuning 
as 5 — /S. — El. In the first step, we nucleate the dipolar 
crystal using a 7r-pulse far detuned from Ei , which corre- 
sponds to 5/^ = VN , where N is the number of sites of 
the ordered structure that one wishes to prepare. Then, 
each of the up spins separated by ~ A^ sites serves as a 
center around which the ordered phase is to be grown. 

Next, we make use of the dipole-dipole interactions be- 
tween the molecules to deterministically grow the com- 
mensurate dipolar crystal with filling 1/n near the nu- 
cleation center. In order to resonantly flip a spin located 
n sites away from the nucleation center we use a tt pulse 
with detuning 5 = V/n^, while frequency flipping the 
third spin in the chain is shifted by (5 = V / v? -\-V / {2nY — 
9F/(8n'^). Such detuned pulses will only affect these two 
spins, flipping any other spin in the system will be highly 
off- resonant. The size of the ordered structure that one 
can deterministically prepare with tt pulses, Axq, is lim- 
ited by the Rabi frequency £7, which must be smaller than 
the change of the detuning 5 from one pulse to another. 
Assuming il ^ 27r x 50 Hz and LiCs molecules, whose 
interaction is given by V/h — 2ti x 40.7 KHz, one can 
deterministically prepare the states with Axq — 8a for 
the n = 2 (AFM) phase and Axq — 9a for the n = 3 
phase. 
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FIG. 3. Pulse sequence used to non-adiabatically prepare 
an ordered phase with a filling 1/n: (i) a far detuned tt pulse 
(blue) nucleates the phase by flipping one spin per A*' lattice 
sites; (ii) a few resonant 7r-pulses (green) deterministically cre- 
ate a few lattice cells of the phase near the nucleation center; 
(iii) continuous driving (red) propagates the crystal bound- 
ary. The required values for the renormalized detunings of 
the pulses, 5 — A — Ei, are shown. The lattice plots exem- 
phfy the 71 = 2 (AFM) case. 



After a few pulses the value of detuning required to 
flip the next spin in the ordered chain approaches the 
value oi S ~ ({3)V/{na)^, and one can use continuous 
driving to grow the ordered structure further. Then, the 
dynamics of the crystal growth is equivalent to a single 
particle hopping on a semi-infinite two-dimensional lat- 
tice {i,j} given by the domain size, Axi > Axq, and 
the domain center of mass position, Xi, as given by the 
effective Hamiltonian, 
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which is schematically shown in Fig. H (a) . Numerically 
solving the time-dependent Schrodinger equation by ex- 
act diagonalization, we study the dynamics of the crys- 
talline domain. As shown in Fig. H (b), the domain size 
increases linearly with time, illustrating that the growth 
process is highly efficient. 

In order to prepare commensurate phases of higher or- 
der, such as 2/5 or 2/7, one can replace the continuous 
driving with a series of pulses, {7ri,7r2,7ri,7r2 . . . }, with 
detunings 5i and 5^ ■ In the case of the 2/5 phase, whose 
unit cell is given by tiiti^ the required detunings are 
5i ~ 0.05 V and 82 ~ 0.14 V . The frequency resolution 
needs to fulfill the condition, Q, < 82 — 5i, which can be 
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FIG. 4. Top: the growth of the crystalline phase can be 
visualized as a particle hopping on a 2D lattice, with Ax the 
domain size, and x the domain center-of-mass position, as 
given by the effective Hamiltonian, eq. ([3|. The example is 
given for Aa;o ~ 3na. Bottom: the dynamics of the growth 
of the ordered phase during the continuous driving. Colors 
indicate the populations corresponding to different domain 
sizes, the inset (same axis labels as in the main figure) shows 
the root-mean-square of the size, which is linearly increasing 
with time. 



easily satisfied for LiCs, requiring 5^ < 4 KHz. 

The leading source of errors in the process of grow- 
ing the dipolar crystal will be due to defects created 
by flipping a spin n -\- 1 sites (instead of n sites) away 
from the domain boundary. During the continuous driv- 
ing, these states are detuned by A„+i = V{1 + [((3) + 
■!/'^^^(l/n)/2]/n'^), where ^^"■\z) is the polygamma func- 
tion [5S]. The number of spins that can be flipped be- 
fore creating a defect can be estimated by the condition 
N^'^/A'^^i sa 1, which effectively hmits the Rabi fre- 
quency. Since the number of spins one can flip within 
the lifetime t of the system is given by A^ ~ fir, we 
can calculate the average size of the domains to be 
N — (A„+it)^/^, which scales more favorably with time 
than the limit set by the Kibble-Zurek prediction. As- 
suming V/h = 27r X 41 kHz and the lifetime of the system 
of T = 1 s, we find that it is possible to prepare domains of 
a n = 2 crystal consisting of up to nN — 1700 spins, while 
for the n = 3 crystal one can reach the size of nN = 1000 
spins. The value of N determines the correlation length 



of the defects, i.e. the effective temperature of the sys- 
tem, Teff ~ A„+i/7V. Here, we find T^s = 200 pK for 
LiCs molecules and N = 1000. 

We note that in principle the performance of this pro- 
cedure can be improved even further, by using a compos- 
ite pulse sequence that dynamically decouples the pro- 
cesses leading to the creation of defects [57]. Likewise, 
employing optical superlattice techniques [21] could also 
be used to improve the addressing of individual transi- 
tions. 

Finally, we would like to emphasize that the dipolar 
crystals prepared after switching off the microwave driv- 
ing correspond to the ground state of the Hamiltonian (pi) 
in the limit oi h^ = 0, and it is possible to explore the 
region of finite hx by adiabatically turning the microwave 
driving field on again. Thus, it becomes possible to ex- 
plore static and dynamic properties of the commensurate 
phases. In particular, the dynamics of fractionalized low- 
energy excitations can be studied, which are eventually 
responsible for the transition to the fioating solid phase 

m- 

In summary, we have proposed a novel method to cre- 
ate strongly interacting many-body states in the dipole 
blockade regime of microwave transitions between rota- 
tional states of polar molecules. Our approach relies 
on the non-adiabatic driving of the microwave transi- 
tions and allows for the possibility to construct larger 
domains of dipolar crystals than in an adiabatic sce- 
nario. The scheme presented here is general and can be 
implemented with any two-state system possessing long- 
range Ising interactions, such as different spin states of 
^S or ^H molecules, trapped ions [29], Rydberg atoms, or 
Nitrogen- Vacancy centers in diamond |23| . Furthermore, 
our proposed method is not limited to dipolar crystals 
and could be applied to other systems where the initial 
state can be efficiently coupled to the many-body state 
of interest by non-adiabatic driving. 
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